Unconventional Fermi surface in two-dimensional systems of Dirac fermions 
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At the low energy regime, the decay rate of two-dimensional massless Dirac fermions due to 
interactions can be written as ImE(aj) oc \uj\^ at zero temperature. We find that the fermion system 
has: I) no sharp Fermi surface and no well-defined quasiparticle peak for < a; < 5; II) a sharp 
Fermi surface but no well-defined quasiparticle peak for \ < x <1; III) both sharp Fermi surface and 
well-defined quasiparticle peak for x > 1. In the presence of long-range gauge/Coulomb interaction 
or certain massless boson mode, the system exhibits unusual behavior belonging to class II. 



o 

(N 



I 

o 
o 



> 

o 

00 

o 
o 



PACS numbers: 71.10.Ca, 71.10.Hf, 71.10.Pm 

Landau's Fermi liquid (FL) theory is the standard 
model of quantum many-body physics During the 
past five decades, this theory has been applied success- 
fully to a wide range of condensed matter systems. In 
any normal FL, there are always sharp Fermi surface 
and well-defined quasiparticle peak at the Fermi en- 
ergy. Generally, the fermion decay rate can be written as 
ImE(aj) cx \uj\^ for small energy uj at zero temperature. 
The quasiparticles are well-defined only when the decay 
rate vanishes faster than uj does upon approaching the 
Fermi surface. Hence, an interacting fermion system is 
a normal FL for x > 1 and a non-FL for x < 1. One 
powerful quantity to characterize the FL is the renor- 
malization factor Z, also called the quasiparticle residue, 
which is defined by the real part of the retarded fermion 
self-energy ReS(w) as 
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It is easy to check that Z = for < a; < 1 but < Z < 1 
for X > 1. Therefore, the residue Z vanishes in a normal 
FL and has a finite value in a non-FL. Physically, Z mea- 
sures the discontinuity of fermion momentum occupation 
number n(k) at Fermi energy, which signals the presence 
of a sharp Fermi surface. 

According to the conventional wisdom, the finiteness 
of residue Z is sufficient to guarantee the presence of a 
sharp Fermi surface and well-defined quasiparticle peak. 
However, the opposite reasoning may not be correct. In 
particular, the presence of sharp Fermi surface does not 
necessarily require that Z is finite. Indeed, there is an 
interesting possibility that a system may have a sharp 
Fermi surface even though Z = 0. Recently, Senthil [2[ 
studied the unusual properties of the quantum critical 
point of a Mott insulator to metal phase transition. In 
his scenario, although Z = at the critical point, there 
can be a sharp Fermi surface, dubbed as critical Fermi 
surface, that is characterized by a kink singularity, i.e., 
the discontinuity in the derivative of the fermion occu- 
pation number 7i(k) at the Fermi energy On the 
experimental side, it has long been known that there is 
Fermi surface (or Fermi arc) in the normal state of un- 
derdoped high temperature superconductors desgite the 
absence of well-defined Landau quasiparticles 



In this paper, we propose that this possibility can be 
realized in some planar interacting systems composed 
of massless Dirac fermions. In recent years, the Dirac 
fermions with linear dispersion have been investigated 
extensively because they are the low-energy excitations 
of many strongly correlated systems, including d-wave 
high temperature superconductor quantum Hall sys- 
tem Q , graphene , and iron-based superconductor Q . 
In different physical systems, these fermions may experi- 
ence gauge, Coulomb, or other kinds of interactions. We 
show that the Dirac fermion system has a sharply defined 
Fermi surface but no quasiparticle peak when the expo- 
nent X of decay rate satisfies ^ < a; < 1. In contrast to 
the scenario of Senthil, the Fermi surface in our case is 
defined by the divergence of the derivative of n(k) at the 
Fermi energy, rather than a kink singularity. When the 
Dirac fermions couple to massless boson modes, such as 
gauge field. Coulomb potential, or certain order parame- 
ter fluctuation, they exhibit such unconventional behav- 
ior. For < X < ^, the fermion system has no sharp 
Fermi surface and no quasiparticle peak, which loses any 
similarity to a free fermion gas. For a; > 1, there are 
both Fermi surface and quasiparticle peak, so the system 
is just a normal FL. 

We begin with the following free action 
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where spinor field 4* describes massless Dirac fermions 
and wf is the fermion velocity. The coupling of Dirac 
fermions with various singular boson modes can result in 
unusual behaviors. 

We first consider the interaction of Dirac fermions 
with an abelian gauge field, which corresponds to three- 
dimensional quantum electrodynamics (QED3). After 
proper modifications, this model can be used to study 
the low-energy properties of d-wave high temperature su- 
perconductor and spin liquid state 11| . Recently, 
we studied the decay rate of massless Dirac fermions due 
to gauge interaction. At zero chemical potential ^ = 0, 
the fermion decay rate is always divergent within pertur- 
bation expansion method. We calculated the decay rate 
by the self-consistent Eliashberg equation approach and 
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found [T2| that it is of the form ImT,{uj) oc at zero 

temperature. At finite fx, the longitudinal component of 
gauge field becomes massive due to static screening, but 
the transverse component remains massless. In this case, 
the decay rate is free of divergence and depends on energy 
as Iml](cj) oc at zero temperature Obvi- 

ously, the chemical potential only changes the coefficient 
of fcrmion decay rate. When applied to cupratc super- 
conductor, the gauge field may couple to an additional 
scalar field In the superconducting ground state, the 
gauge field acquires a finite mass via the Anderson-Higgs 
mechanism. It is easy to obtain that, Iml](aj) oc 
which has the same energy dependence as the decay rate 
yielded by contact four-fermion interaction [l5| . 

In the context of graphene, the massless Dirac fermions 
experience the Coulomb interaction, which is unscreened 
due to the vanishing density of states at Fermi energy Q. 
From previously analysis [T6j , we know that the decay 
rate caused by Coulomb interaction is Iml](cj) oc |a;|, 
which is clearly marginal FL behavior [l^l . 

In the vicinity of continuous quantum phase transition, 
the massless Dirac fermions interact with the strong fiuc- 
tuation of order parameter [isj . Such interaction can be 
described by a Yukawa coupling between a spinor field 
and a scalar field. In particular, at the critical point 
between two su per condcuting phases in d-wave cuprate 
superconductor [19[ , the fermion decay rate due to mass- 
less order parameter fluctuation is Iml](a;) oc 

In general, after including the self-energy correction 
due to various interactions, the full retarded Green func- 
tion of Dirac fermion is 



G(w,k) 
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where = — ii with /i being the chemical potential 
and Ek — fF|k| energy of the Dirac fermion. To sim- 
plify the problem, here we assume that the self-energy 
depends only on energy lj. The spectral function, de- 
fined as ^(cj,k) = — ilmG(w,k), has the form 



A{LU,k)= — 



Iml](w) 



1 



(4) 



The fermion momentum occupation number is given by 



n(k) 



dujA (w, k) . 



(5) 



In order to make a general analysis, we write the zero- 
temperature decay rate of Dirac fermions as 



Iml](cj) =C|wr, 



(6) 



where x is a tuning parameter and C a negative constant. 
The real part of self-energy can be obtained from the 
Kramers-Kronig relation 



Rel](w) = \P j 
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We first consider the region < x < 1. In such region, 
the real self-energy function is found to be 



Rel](w) = Csgn(w)|w|"=/(a;), 



(8) 



where I{x) is 
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When < X < 1, the function I{x) has a finite mag- 
nitude and hence the real part of retarded self-energy 
ReS(a;) has the same energy dependence as ImE(w). In 
particular, we have /(1/2) = 1 and /(2/3) = y/S. 

When X > 1, there will be divergence when Eq.Q 
is used to calculate ReE(aj). To get a finite result, we 
should introduce a cutoff energy LUf and write 



ReE (oj) = -P 

TT 

For a; = 1, we have 
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ReE(w) = Csgn(cj)|cj|- 

TT 
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which in the low-energy regime reduces to 
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This corresponds to the marginal FL behavior. For x > 
1, the real part of self-energy at low energy is 



ReE(w) w G ui. 



(12) 
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FIG. 1: Spectral function, (a) < x < i; (b) | < a; < 1; (c) 
1 < s < 2; (d) X > 2. 

Using these self-energy functions, it is easy to obtain 
that Z = OforO<x<l and < Z < 1 for a; > 1. From 
the spectral functions at various x shown in FigHJ we see 
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that the well-defined quasiparticle peak exists only for 
X > 1 where Z ^ 0. The residue Z also represents the 
discontinuity of momentum occupation number n(k) at 
|k| = hp. When < x < 1, there is no discontinuity in 
n(k); but when x > 1, there is a discontinuity in 7i(k). In 
conventional many-body theory, Z can uniquely differen- 
tiate a non-FL from the normal FL. In strongly correlated 
systems, however. Z itself can no longer fully character- 
ize the non-FL. In order to see how this happens, we will 
calculate n(k) and its derivative for all values of x. 
For < X < 1, the momentum occupation number is 



.(k) 



At Fermi energy, its derivative with respect to |k| is 



n'(k)|, 



oc 



duj- 



As w — 0, the integrand approaches — . For < a; < 
^, it is easy to obtain 

"'(k)lfc^_o = finite, n'(k)|j^^^o = finite. (13) 

For 5 < x < 1, '^'(k) diverges as: 

n'(k)|^^_o ^ -oo, 7i'(k)|j^^+o ^ -oo, (14) 

which indicates that the momentum occupation number 
drops dramatically at fcj?. 

Similarly, for x = 1, we can obtain the following deriva- 
tive of occupation number 



dhJ- 



-C^wln 



cj - C^uj In 



As w — > 0, the integrand approaches — 
-oo, n'(k)|j^^_^o 



^{C\uj\r 

73-, thus 



2h is similar 
For X > 1, we have 



(15) 



which is similar to the case of | < x < 1. 
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0, the integrand approaches — , 



For 1 < 



X < 2, the residue Z 7^ and the derivative of 7i(k) is 
For X > 2, although Z 7^ 0, we know that 



n'(k)|j^^_o = finite, n'(k)| 



feF-l-O 



finite. 



(17) 
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FIG. 2: Fermion momentum occupation number 7i(k). (a) 
< 2: < i; (b) i < 2: < 1; (c) 1 < a; < 2; (d) a; > 2. 
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FIG. 3: Derivative of n(k). (a) < a; < i; (b) i < a: < 1; (c) 
1 < x < 2; (d) a; > 2. 



The dependence of momentum occupation number 
n(k) on X is shown in Figl^l Apparently, as the param- 
eter X grows, the interacting system of massless Dirac 
fermions first develops a sharp Fermi surface at x = i 
and then develops a finite quasiparticle residue Z once x 
exceeds unity. There is a fundamental difference between 
the n(k) functions for < .x < i and for i < x < 1. In 
the former case, the residue Z = 0, so there are no well- 
defined quasiparticles. Moreover, the derivative of ?T.(k) 
is continuous at the Fermi energy kp, so there is also no 
Fermi surface. In the latter case, the residue Z ~ and 
thus there is no well-defined quasiparticle peak. How- 
ever, the derivative of n(k) diverges at kp, therefore n(k) 
drops suddenly as |k| increases across kp, which can be 
identified as the presence of a sharp Fermi surface. To 
see more details of the evolution of n(k) with x, we also 
show the derivative of n(k) in FigO 

In summary, we can divide all the interacting Dirac 
fermion systems into three classes: I) for < x < ^, 
there are no sharp Fermi surface and no well-defined 
quasiparticle peak; II) for ^ < x < 1, there is sharp 
Fermi surface but no well-defined quasiparticle peak; III) 
for X > 1, there are both sharp Fermi surface and well- 
defined quasiparticle peak. 

In QED3, the fermion decay rate has x = ^ at zero /x 
and X = I at finite /i. These states fall into the class 
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FIG. 4: n(k) for different chemical potential fi in QED3. is 
zero in (a) and increases from (b) to (d). 




FIG. 5: Critical Fermi surface of Senthil. (a) n(k); (b) n'(k). 

II. To see the Fermi surface evolution with growing /i, 
we show n(k) in Fig|4l We can see that n(k) does not 
evolve continuously from /i = to the case of finite /i. On 
the contrary, as /x jumps from zero to a small value, the 



curve of n(k) changes dramatically. However, the Fermi 
surface remains sharply defined for arbitrary small value 
of fi, because u does not alter the energy dependence of 
decay rate The marginal FL (a; = 1) caused by the 
long-range Coulomb interaction or certain massless or- 
der parameter fiuctuation also belongs to class II. In the 
presence of short-ranged gauge/Coulomb interaction, the 
decay rate has exponent x > 1 and thus the behavior falls 
into class III. Unfortunately, it is currently unclear what 
kind of interactions can produce behaviors belonging to 
class I, which deserves further research. 

We emphasize that the unconventional Fermi surface 
studied in this paper is different from that of Senthil ^ . 
As shown in FiglSJ the critical Fermi surface proposed by 
Senthil is defined by the discontinuity in the derivative 
of n(k). In our problem, however, there is no such kink 
singularity and the sharp Fermi surface is characterized 
by the divergence of the derivative of n(k). 

Finally, although all the above discussions were made 
in the case of massless Dirac fermions, the unconventional 
Fermi surface of class II should also be realized in a num- 
ber of non-relativistic strongly correlated systems, such 
as non-relativistic fcrmion-gauge system with large Fermi 
surface [20j and quantum critical metals near various 
symmetry-breaking instabilities [2l|, [22|, provided that 
the exponent in fermion decay rate satisfies < x < 1. 
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